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FOREWORD 
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for  Potential  Chemical  Agents  (U).  This  work  was  started  in  March  1969  and  completed  in  October 
1970.  The  data  are  recorded  in  notebook  8463. 

Reproduction  of  this  document  in  whole  or  in  part  is  prohibited  except  with  permission 
of  the  Commanding  Officer,  Edgewood  Arsenal,  ATTN:  SMUEA-TSTI-T,  Edgewood  Arsenal, 
Maryland  21010;  however,  DDC  is  authorized  to  reproduce  the  document  for  United  States 
Government  purposes. 

The  information  in  this  document  has  not  been  cleared  for  release  to  the  general  public. 
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DIGEST 


A  least-squares  technique  has  been  developed  for  determination  of  the  constants  of  the 
Antoine  equation.  The  method  involves  use  of  the  procedure  of  search  by  golden  section.  A 
conversational  FORTRAN  program  is  provided  with  results  of  sample  runs  utilizing  data  from  the 
open  literature.  The  method  gives  accurate  fits,  though  the  constants  calculated  differ  slightly  from 
the  literature  constants  because  the  method  described  is  for  fitting  data  that  have  a  relatively 
constant  percentage  error  rather  than  a  constant  absolute  error.  The  former  is  usually  more 
appropriate  when  values  are  measured  over  a  wide  range  of  temperatures.  Plots  are  provided 
automatically  through  the  use  of  a  teleplotter. 
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CONVERSATIONAL  COMPUTATION  METHOD  FOR  FITTING  THE  ANTOINE  EQUATION 
TO  VAPOR  PRESSURE-TEMPERATURE  DATA 


I.  INTRODUCTION. 

The  vapor  pressure  of  a  pure  compound  is  a  unique  function  of  temperature.  It  is  a 
property  required  for  many  calculations  in  physical  chemistry,  chemical  engineering,  and  numerous 
other  disciplines.  Equations  relating  to  vapor  pressure  and  temperature  are  required  to  summarize 
vapor  pressure  data;  to  permit  the  use  of  the  data  in  computer  programs;  and  for  smoothing, 
interpolation,  extrapolation,  and  determining  heats  of  vaporization.  One  of  the  most  widely  used 
equations  for  fitting  vapor  pressure-temperature  data  is  the  Antoine  equation,1  the  use  of  which  has 
been  reviewed  by  Thomson.2 

In  addition  to  fitting  vapor  pressure  data,  the  Antoine  equation  may  be  used  to  fit 
densities,  fugacities  of  liquids,  kinematic  viscosities,  relative  vapor  pressures,  and  other  types  of 
data.  The  Antoine  equation,  which  has  its  basis  in  the  well-known  Clapeyron  equation  (see  any 
elementary  physical  chemistry  text),  is* 

dP  AH v  (i) 

dT  TAV 

where 

P  =  vapor  pressure 
T  =  absolute  temperature 
A Hv  =  heat  of  vaporization 

A  V  =  change  in  volume  associated  with  vaporization 

If  the  molar  volume  of  the  liquid  is  negligible  compared  with  that  of  the  vapor  and  the 
vapor  is  a  perfect  gas,  then  the  Clausius-Clapeyron  equation  is  applicable. 

d\nP  _ 

cKUT)  '  R  (2) 

This  equation  may  be  integrated,  if  AHV  is  constant,  to  give 

\ogP  =  A  -  B/T  (3) 

where  A  and  B  are  parameters.  The  equation  also  may  be  written  as  follows: 

logP  =  A  -  B/(t  +  273.16)  (4) 


where  t  is  temperature  in  °C. 

The  Antoine  equation  results  from  replacing  the  value  273.16  in  equation  (4)  with  a 
parameter  C 


*A11  terms  are  defined  in  the  glossary  at  the  end  of  the  text. 
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log  P  =  A-  Bl(t  +  Q 


(5) 


This  equation  has  been  found  to  be  adequate,  except  for  data  of  very  high  precision,  over  the  range 
from  the  triple  point  to  reduced  pressures  of  0.75. 2 

The  advantages  of  the  Antoine  equation  include  the  following:  small  number  of  constants 
required,  temperature  may  be  made  the  dependent  variable  by  a  simple  inversion,  and 
extrapolations  are  believed  to  be  more  reliable  than  for  other  equations.3 


The  value  for  C  is  usually  found  to  be  lower  than  273.15  and  has  often  been  taken 
arbitrarily  as  230.4  The  latter  value  has  been  found  to  be  typical  for  organic  compounds  that  are 
liquids  at  room  temperature.  Because  calculating  the  best  value  for  C  is  a  time-consuming  task 
without  a  computer,  the  value  C  =  230  frequently  is  found  in  the  literature. 


Equations  (2)  and  (5)  can  be  combined  to  give  the  heat  of  vaporization 

Afi'=23RB(jc)2 

II.  THEORY. 


(6) 


The  Antoine  equation  [equation  (5)]  may  be  written  in  the  form  of  a  linear  equation 

Y  =  A+BX  (7) 


by  use  of  the  transformations 


Y=\ogP  (8) 

X  =  -\/(C+t)  (9) 

A  set  of  Y’s  and  Y’s  may  be  fitted  by  ordinary  linear  regression  analysis,5  which  results  in  the 
following  equations: 


a  =  (  E  r,-n  £  x\ 

\  i=l, h  f=H,n  / 


In 


(10) 


B 


■£z’-CS4i 


do 
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where  Yt  and  Xt  are  the  ith  value  of  Y  and  X,  respectively,  and  n  is  the  number  of  points  being 
fitted.  Given  a  set  of  vapor  pressures  with  the  corresponding  temperatures,  any  value  of  C,  4 ,  and  B 
may  be  found. 

The  method  of  standard  error  of  estimate  involves  minimizing  the  following  function, 


£[ 

i=l,  n 


Yj  (experimental)  -  (calculated  from  eq 


6)3  2 


which  will  be  referred  to  as  S.  In  this  case,  S  may  be  considered  a  function  of  C  alone  because  both 
4  and  B  are  known  functions  of  C  for  a  given  set  of  experimental  data.  Because  S  may  not  be 
minimized  analytically  by  any  simple  known  method,  values  of  S  are  calculated  for  varying  values 
of  C  until  a  minimum  value  of  S  is  found. 

The  following  equation  results  from  the  above  discussion: 

S  =  2  IP, -4 -£*,-]  2  (12) 

1=1,13 


Expansion  yields 


5  = 


2  (A2  +  2 ABXt  +  B2Xt2  +  Y,-2  -  2A  Y{  -  2 BX,Yt) 

i=l,  n 


(13) 


Summing  terms  individually  gives 


S  =  nA2  +  2 AB  La  Xt  +  B2 

3  “1,13  1=1,13 


2  Xf2  +  2  r,-2  -  24  £  Yi  -  2 B  £ 


/=1,  n 


Yt2  -24  La  Yr2B  La  X,  Yi  (14) 

1=1,13  1=1,13 


which  is  used  in  the  computer  program  to  determine  S. 

The  search  method  used  to  find  the  best  value  of  C  involves  division  of  a  line  segment  by 
a  method  developed  by  Euclid.  The  method,  called  the  search  by  golden  section  procedure  and 
described  in  detail  by  Wilde,6  had  great  mystical  significance  among  scholars  of  the  Renaissance. 
This  search  method  was  selected  because  it  does  not  require  advance  knowledge  of  how  many 
calculation  steps  are  required.  It  involves  division  of  a  line  segment  into  two  unequal  parts  so  that 
the  ratio  of  the  whole  to  the  larger  part  is  equal  to  the  ratio  of  the  larger  to  the  smaller  segment. 
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In  mathematical  terms  this  definition  leads  to  the  following  two  equations: 


First, 


Lj  Lj+ 1 


+  Lj+2 


where 


Lj  =  length  of  a  line 

I  Lj+  j  ~  1  2  =[length  of  segments  into  which  the  line  is  divided 
Second, 


^7+ 1  ^/+2 

It  may  readily  be  shown  from  the  above  equations  that  Z  must  be  equal  to  1.6180340  and  that 

Lj+ 1  =4^-=  0.61803399  Lj 

After  m  divisions  of  the  line  (which  is  equivalent  to  dividing  the  original  line,  L0 ,  m  times) 


rn  jjn 

As  an  example,  after  24  divisions  of  the  line,  Lm  is  equal  to  0.00000965  times  the  length  of  the 
original  line,  La . 

The  standard  error  of  estimate  provided  by  the  program  is  calculated  from  the  following 

equation: 

<i5) 

It  is  emphasized  that  the  values  for  the  logarithm  of  vapor  pressure  are  used  for  the 
dependent  variable  in  the  least-squares  fit;  the  calculated  logarithms  of  vapor  pressures  are 
subtracted  from  the  logarithms  of  the  measured  pressures  to  calculate  the  standard  deviation.  Using 
the  data  in  this  manner  prevents  excessive  weighting  of  larger  pressures  in  the  analysis.  Generally, 
the  method  tends  to  weight  values  as  a  fixed  percentage  of  the  measured  value.  If  the  absolute  error 
of  the  measurement  is  constant,  regardless  of  the  magnitude  of  the  values  measured,  then  the 
method  developed  in  this  report  probably  will  not  provide  the  best  possible  fit  to  the  data. 

III.  DESCRIPTION  OF  COMPUTER  PROGRAM. 

The  program  is  quite  self-explanatory;  the  first  question  the  program  asks  is  how  many 
pairs  of  temperature-pressure  values  are  to  be  entered.  This  is  answered  by  typing  the  number  of 
data  pairs  plus  a  carriage  return.  Next,  the  points  are  entered;  as  many  as  100  points  may  be  used, 
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but  only  70  characters  per  line  may  be  entered  and  recorded  by  the  teletype.  There  must  be  a 
carriage  return  at  the  end  of  each  line.  If  the  points  are  to  be  entered  by  paper  tape,  at  the  end  of 
each  line  there  must  be  a  carriage  return,  X-off,  and  rubout. 

The  computer  queries  the  operator  regarding  the  need  for  a  data  table,  after  which  the 
equation  constants  are  printed  along  with  the  calculated  standard  deviation  and  an  extrapolated 
boiling  point.*  This  information  is  followed  by  a  printout  of  the  experimental  and  calculated  vapor 
pressure  values. 

An  option  then  is  provided  for  generating  a  table  of  vapor  pressures,  volatilities,  and  heats 
of  vaporization  over  a  temperature  range  and  in  temperature  increments  selected  by  the  operator. 
This  information  is  printed  after  the  molecular  weight  of  the  compound  is  provided. 

The  program  automatically  loads  the  experimental  and  calculated  values  of  vapor  pressure 
into  a  file  named  EZ  for  plotting  on  a  Typagraph  Corporation  teleplotter.7 

Generally,  when  vapor  pressures  have  been  measured  only  over  a  narrow  range  of 
temperatures  or  the  data  are  not  very  accurate,  it  is  better  to  estimate  a  value  for  C  and  then  do  a 
standard  linear  least-squares  fit  of  the  data.  When  C  is  calculated  and  found  to  be  greatly  different 
from  similar  compounds  of  comparable  molecular  weight,  one  should  be  suspicious  of  the  data. 
Thomson2  has  discussed  these  matters  in  detail. 

Appendix  A  defines  the  items  used  in  the  computer  program,  which  is  given  in  appendix 
B. 

IV.  SAMPLE  CALCULATIONS. 

Appendix  C  provides  two  sample  runs  using  the  reported  vapor  pressure  data  of  Kremme 
and  Kreps.8  The  first  sample  run  is  for  1  -hexadecanol.  The  vapor  pressures  for  this  compound  were 
plotted  automatically  and  are  shown  in  the  figure.  Although  Kremme  and  Kreps  used  a  different 
method  of  determining  the  Antoine  constants,  they  are  in  good  agreement  with  the  values 
calculated,  as  shown  in  table  I.  The  parameters  differ  by  1%,  and  the  average  and  maximum  errors 
of  estimate  are  comparable. 

The  second  sample  run  is  for  1  -tetradecanol.  This  sample  run  demonstrates  a  situation  in 
which  C  is  found  to  be  outside  the  range  from  100  to  273.16.  In  this  case,  the  constants  differ 
significantly  from  those  determined  by  Kremme  and  Kreps.8 

It  was  found  that  the  Antoine  constants  of  Kremme  and  Kreps  led  to  different  errors  of 
the  estimate  than  were  reported.  Recalculation  of  the  errors  with  their  data  and  constants  shows  a 
maximum  error  of  1 1.4%  instead  of  the  reported  value  of  1.665%.  Table  II  provides  the  results  of 
the  recalculation. 

This  computer  program  was  run  also  with  the  vapor  pressure  data  of  isobutyl  borate 
reported  by  Christopher  and  Shilman.9  The  calculated  boiling  point  of  210.4°C  agreed  exactly  with 
their  value.  They  calculated  a  heat  of  vaporization  of  12.345  kcal/mole,  which  compared  favorably 
with  the  value  of  12.33  obtained  in  this  study. 


•Caution  should  be  exercised  in  using  extrapolated  vapor  pressure  values.  Depending  on  the  change  in  heat  of  vapor¬ 
ization  with  temperature  and  the  extent  of  extrapolation,  significant  errors  may  result. 
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+1*740  +1*880  +2*020  +2*160 
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+2*300 
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+  13 
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*  13 

-?  CALCULATED 

Figure.  Typagraph  Plot  of  Vapor  Pressure  Versus  Reciprocal 

Absolute  Temperature  for  1-Hexadecanol 

Table  I.  Comparison  of  Data  Fits  by  Kremme  and  Kreps*  With  Sample  Runs 


Compound 

Antoine  constants 

Error  of  estimate 

i  Calculation 

A 

B 

C 

Avg 

Max 

'  method 

1 1-Hexadecanol 

7.04183 

1880.126 

127.123 

e, 

1.180 

b 

3.569 

Kremme  &  Kreps. 

7.06077 

1893.76 

128.406 

1.189 

3.536 

This  report. 

1-Tetradecanol  , 

6.48407 

1412.907 

95.368 

0.793 

1.665 

Kremme  &  Kreps. 

6.21962 

1244.90 

75.600 

1.6 

2.072 

11.4 

5.775 

Recalculated. 

This  report. 

♦Kremme,  H.  R.,  and  Kreps,  S.  I.  J.  Chem.  Eng.  Data  14,  98  (1969). 


Table  II.  Recalculated  1-Tetradecanol  Data  I 


Temp 

Vapor  pressure 

Difference 

Exper 

Calcd 

°C 

Torr 

% 

151.6 

5.2 

5.8 

11.4 

163.0 

10.5 

10.4 

-1.3 

171.0 

15.3 

15.1 

-1.1 

177.0 

19.9 

19.8 

-0.5 

188.4 

32.1 

32.0 

-0.4 

199.1 

48.4 

48.5 

0.2 

213.3 

80.1 

80.7 

0.7 

225.7 

120.0 

121.2 

1.0 

243.1 

202.8 

204.0 

0.6 

257.8 

302.6 

304.4 

0.6 

277.9 

501.2 

499.8 

-0.3 

295.9 

755.2 

_ i 

746.4 

-1.2 

The  program,  with  slight  modification,  may  be  used  to  fit  other  physicochemical  data; 
i.e.,  densities,  fugacities  of  liquids,  kinematic  viscosities,  etc. 

V.  CONCLUSIONS. 

A  conversational  FORTRAN  computer  program  is  described  that  calculates  the  constants 
required  for  solution  of  the  Antoine  equation  for  vapor  pressure-temperature  data.  Convenient 
printouts  of  experimental  data  used,  calculated  vapor  pressure  data,  volatility  data,  heats  of 
vaporization,  and  an  estimated  normal  boiling  point  are  provided.  In  addition,  the  standard  error  of 
estimate  is  calculated. 
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GLOSSARY 

A,B,C  Parameters  of  the  Antoine  equation 
kHv  Heat  of  vaporization 

Lj  Line  segments  where  /  is  number  of  divisions  by  golden  section  method 
n  Number  of  data  points 
P  Vapor  pressure 

A  Ideal  gas  law  constant 

S  Sum  of  the  squares  of  differences  between  calculated  and  measured  values  of  log  P 
SD  Standard  error  of  estimate 
t  Temperature  in  °C 

T  Temperature  in  °K 

A  V  Change  in  volume  asspciated  with  vaporization 
X  Reciprocal  of  -  (C  +  f ) 

Xf  The  ith  value  of  X 

Y  Logarithm  of  P 

Yj  The  ith  value  of  Y 

Z  Ratio  of  successive  lengths  by  golden  section 
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A,  B 

BP 

B9 

C 

Cl 

C2 

D1 

G 

GOLD 

H 

I 

IC 

INT 

II,  J,K,K1 

L 

LI 

N1 

N9 

PCD 

PRC 

RDT 

SD 

SI 


APPENDIX  A 

DEFINITIONS  OF  TERMS  IN  PROGRAM 
Parameters  in  the  Antoine  equation 
Boiling  point  calculated  from  Antoine  equation 

Maximum  allowable  range  of  uncertainty  associated  with  the  value  of  parameter 
C 

Parameter  in  the  Antoine  equation 
Log  P  calculated  from  Antoine  equation 
P  calculated  from  Antoine  equation 
Controls  printing  of  optional  data  table 
Linear  array  of  measured  pressure  values 
Line  segment  from  golden  section  method 
Maximum  temperature  used  in  optional  output  table 
subscript 

Subscript  to  be  used  in  sub-routine  names  COMP  to  select  C 

Interval  of  temperature  for  calculating  table 

Subscripts 

Subscript  or  lowest  temperature  in  optional  output  table 
Subscript 

Number  of  data  points 

Controls  printing  of  diagnostic  message 

Percent  difference  between  measured  and  calculated  pressures 

Pressure  calculated  from  Antoine  equation 

Reciprocal  of  temperature  in  °K 

Standard  error  of  estimate 

Sum  of  Xj’s 
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52  Sum  of  X42  ’s 

53  Sum  of(Xj  *  Yj)’s 

54  SumofYj’s 

55  Sum  of  Yj2’s 

T  Linear  array  of  experimental  temperatures  in  °C 

T1  Temperature  used  for  calculating  optional  output  table 

V  Linear  array  of  successive  values  of  S 

W  Linear  array  of  trial  values  of  parameter  C 

W1  Lowest  value  of  C  used;  normally  100 

W2  Highest  value  of  C  used;  normally  273 

X  Reciprocal  of  -  (C  +  t) 

Y  Linear  array  of  common  logarithms  of  measured  pressure  values 
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APPENDIX  B 


PROGRAM  LISTING 

The  language  used  is  CALL/360:  FORTRAN**.  CALL/360  is  a  system  designed  to 
provide  the  enormous  efficiency  of  fast  turnaround  conversational  computing  with  immediate 
access  to  the  full  capabilities  of  a  large  computer,  IBM  360/65  in  this  case,  at  a  conveniently  located 
remote  terminal.  Because  a  few  of  the  features  of  this  subset  of  FORTRAN  may  not  be  familiar  to 
some  readers,  they  are  discussed  below. 

First,  the  REAL*8  specification  statement  establishes  variables  as  real  constants  of  eight 
locations  or  approximately  16.8  decimal  digits  of  precision.  The  rest  of  the  real  floating  point 
variables  are  about  7.2  decimal  digits.  Second,  the  percentage  symbol  at  the  end  of  a  line  replaces 
the  continuation  number  commonly  used  with  cards.  Third,  the  “*”  used  in  the  READ  commands 
in  place  of  a  format  statement  number  indicates  that  the  format  to  be  used  for  data  input  is  list 
directed.  In  this  format,  data  items  are  separated  by  either  a  comma  or  a  blank  space. 


**CALL/360  FORTRAN  Language  Reference  Manual.  G320-1054-0.  International  Business  Machines  Corporation, 
New  York,  New  York.  1970. 


JCP-I,  DPG 


1  00  "TI1LH  AN  10  INE 

105  “PROGRAM  TO  FIND  CONSTANTS  FOR  THE  ANTOINE  EQUATION 
110  COMMON  Y( 1 00) #  T( 100)* S4# S5# A*  F*  W(4)#V(4)*N 1  #  I C 

1  15  DIMENSION  G< 100) 

180  HEAL*8  Y*T#S4*S5#A*B*  AFL 

125  CALL  OPEN<  1*  *  EZ  '*  'OUTPUT  ) 

1 30  “INPUT  SECT  ION 
135  WRI  TE(  6#  1  00) 

140  100  FORMAT  </*  UNITS*  TEMPERATURE  CENTIGRADE*/* 

145  *  PRESSURE  TORR’/) 

1 50  1  WRITE<6* 10) 

155  10  FORMAT  <*  HOW  MANY  FOINTS?') 

160  READ(  5#  *  )N  1 

165  S4»0 

170  XXX-1E37 

175  XXL-9E37 

180  S5*0 

185  WRI  T  E<  6#  1  1 ) 

190  11  FORMAT  (*  DATA*  T1*P1*T2*P2# . *) 

195  READ(  5**)<T(I)#G(I)#I»1#N1) 

2  00  “CALCULATE  LOGS  OF  VAPOR  PRESSURES#  SUM  OF  LOGS*  AND  SUM 
2  05  “OF  LOGS  SQUARED 

2  10  DO  20  I»  1*  N 1 
2  15  RTK*1. 0/(27  3.  16+7(1)) 

2  20  WRI  T  E(  1)  RTK*  G(  I  ) 

225  Y( I )*AL0G10(G( I ) ) 

2  30  S4*S4+Y (  I  ) 

2  35  S5*S5+Y( I )*Y( I ) 

240  20  CONTINUE 
2  45  WRITE(l)  XXX*XXX 

2  50  140  WRI 1E(6*  12) 

2  55  12  FORMAT(*  DATA  TABLET  (  1-YESl  0*N0  )  *  ) 

2  60  READ<5#*>  D1 

265  41  FORM  AT  (9X*  TEMPERA!  URE*6X*  PRESSURE*/ > 

270  I F  (  D1  .LT.  0.5  )  GO  TO  156 

275  WRI  1  E(  6*  4  1 ) 

280  DO  40  I«l#  N 1 

285  WRI  T  E(  6#  87  )  1(I)#G(I) 

290  87  FORMAT  (  3X*  1  PE  16.  4*  1  PE  16*  3) 

295  40  CONTINUE 

3  00  “DEFINE  RANGE  OF  SEARCH 
3  05  156  W 1-1 00 

3  10  W2-273*  15 

3  15  “  BEGIN  SEARCH 
320  160  CONTINUE 

3  25  “SET  LIMIT  OF  UNCERTAINTY  (B9)  ALLOWABLE  FOR  C 
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3  30  E9*0*  005 

335  1-1 

3  40  d-2 

345  K* 3 

350  L-4 

3  55  "THE  FOUR  POINTS  INVOLVED  IN  THE  SEARCH  ARE  STORED  IN 
3  60  "  V<N>«  N-l#4  IN  AN  ORDEP  DETERMINED  BY  PROGRESS 

365  "  OF  SEARCH 
370  WCD-V1 

375  IC-I 

3  80  CALL  COMP 

385  WCL>* W2 

390  I  C*L 

395  CALL  COMP 

400  184  GOLD*fcCL)-KI> 

4  05  GOLD* GOLD*  •  6  1803399 

4  10  IFC GOLD.LT*  E9>  GO  TO  258 

4  15  VCKJ-WCI  )+GOLD 

4  20  I  C*K 

4  25  CALL  COMP 

4  30  "IF  STATEMENTS  INVOLVING  V(N)  DETERMINE  WHICH  LINE  SEGMENTS 
4  35  "TO 

440  I F  C  V(  I  )  •  GT#  V(L  )  )  GO  TO  202 

445  L-l 

450  1*4 

455  202  IFC  VCIO.LT.VCL))  GO  TO  214 

460  L*4 

465  WCD-VCK) 

470  VC  I  )  *VCK) 

475  1*1 

4  80  GO  TO  184 

485  214  L*4 

490  1-1 

495  218  GOLD-GOLD*.  6  1803399 

5  00  I F (  G0LD.LT.B9)  GO  TO  258 

505  WC  JJ-VCI  )*GOLD 

5  10  IC-J 

5  15  226  CALL  COMP 

520  IFC  V(J).LT.VCK))  GO  TO  248 

525  I  1- 1 

530  I-J 

5  35  J-K 

5  40  K* I  1 

5  45  IFC  G0LD.LT.E9)  GO  TO  258 

550  WCK)-WCI  )+GOLD 

5  55  GOLD- GOLD*.  6 1803399 

560  I C-K 

5  65  GO  TO  226 
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5  70  2 48  Ll-L 

575  L=K 

580  K-J 

585  J-Ll 

590  GO  TO  218 

595  "DETERMINE  WHETHER  C  IS  IN  THE  INITIAL  RANGE 
600  258  IFtWl.EG.W2)  GO  TO  301 

6  05  IC«J 

6  10  I F  t  W(I).GT.Wl)  GO  TO  266 

6  15  IC-I 

6  20  GO  TO  356 

625  266  I  F  t  VtL>.LT.W2>  GO  TO  300 

6  30  IC«L 

6  35  GO  TO  362 

640  301  IC-I 

6  45  300  CONTINUE 

650  C-W(IC) 

6  55  WRI  TE(  6#  29  ) 

6  60  29  FORMATt  '  0  LOG  P- A-B/t  C+T)  '  ) 

665  WRI  TE(  6*  28  >  A#B#WtIC) 

6  70  28  FORMATt  '  0  A*  5X#  1PE 1 2.  5/ 3X  ' B  '  5X#  1PE  12.  5/3X  *C  *5X#  1PE  1 2.  5) 
675  SD-SQETtVt  IC)/tNl- 1)) 

680  WRI  TEt  6#  27  )  SD 

6  85  2?  FORMATt  *  0  STANDARD  DEVI  ATI ON ' 3X#  1PF12.  5) 

690  BP«-B/(2.88081-A)-C 

695  WRI  TEC  6#  26)BP 

700  26  FORMATt  *0  BOILING  FOINT  '  4X#  1PE12.5) 

7  05  WRI  TL(  6#  38  > 

710  38  FORMATt ' OTEMPERATURE ' 2 IX 'PRESSURE ' 18X' PERCENT '/* 

7  15  22X 'EXPERIMENTAL  '5X 'CALCULATED' 1  OX ' DI FFERENCE '  ) 

720  DO  111  I-1#N1 

7  25  PRC- A-E/  t C+Tt I  )  ) 

7  30  PRC«10.**PRC 

7  35  PCD- 1  PRC-Gt  I  )  )*  1 00*  0/6t  I  ) 

740  43  FORMAT  t  2X#  F8«  3#  1 2X#  F 1  1  •  5#  5X#  F 1  1 . 5#  6X#  F 1 1 . 5) 

7  45  RTK-  1*  0/  t  273*  16+Tt I ) ) 

7  50  WRI  TEt  1 )  RTK#  PRC 

7  55  111  WRI  TEt  6#  43)  Ttl)#  Gtl>#  PRC#  PCD 
7  60  WRI  TEt  1)  XXL#  XXL 

7  65  CALL  CLOSEtl) 

770  39  FORMATt '0  TABLE?  GIVE!  LOWEST#  HI  GHEST  AND  INCREMENT  OF* 
775  TEMPERATURE'/'  AS  INTEGERS  t  0#  0#  0  MEANS  NO)') 

7  80  WRI  TEt  6#  39) 

785  READt  5#*  )  L#H#  INT 

7  90  I  Ft  L.EG.H)  GO  TO  387 

795  WRI  TEt  6#  101 ) 

8  00  101  FORMATt/'  GIVE  MOLECULAR  WEIGHT') 

8  05  READt  5#  *  >  XMOL 
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8  10  UhIU(6>46) 

8  15  46  FORMAT'  1  EMPERA1 URE  VAPOh  PRESSURE  VOLATILITY* 

8  20  HEAT  Oh  VAP0RIZA1  ION  *  ) 

8  25  1  1=L 

8  30  h=0. 082053 

8  35  C0N= 1.9872E-3+2. 3026 

840  344  Cl=A-b/(C+l  1) 

845  346  IK  1  l.GT.H)  GO  10  388 

8  50  C2=10.  0+  +  C1 

8  55  18=1  1  +  273.  15 

8  60  V0L=C2*XM0L/  Ch+lK+760.  0)*  1. 0E.+  06 

865  HEA1=C0N+E*< IK/C  1  1  +  C) )**2 

870  WhI  TE(  6.  59 )  11.  C2*  V0L*HEA1 

875  59  E0hMAK2X*  E  1  0.  3*  E  1  6.  5*  E  1 6.  3*  F  1 6.  3 ) 

880  1  1=1  1  +  IN  1 

8  85  GO  10  344 

890  356  WhI 1 E( 6*  49 )  W(IC) 

895  49  FORMAT  '0  C  LESS  1HAN  '  1PE1&.5) 

9  00  GO  10  366 

905  362  WRI'1E(  6.  48  )  W(IC) 

9  10  48  FORMAIC'O  C  GREAT  Eh  THAN  '1PE12.5) 

9  15  366  IKN9.G1.0)  GO  10  378 

9  20  N9=  3 

925  WhI  1  E(  6.  47  ) 

930  47  FORMAT '01HE  CONSTANT  C  IN  LOG  P*A+B/!C+1)  IS  ESTIMATED  10  BE’/X 
9  35  •  OUTSIDE  THE  RANGE  100  10  273.  10  CONTINUE*  SPECIFY  THE'/X 

9  40  *  ENDPOINTS  OE  AN01HER  hANGE.  IE  BOTH  ENDPOINTS  ARE  THE*/* 

9  45  '  SAME*  THE  CALCULATION  WILL  BE  DONE  USING  1HAT  VALUE  FOR'/X 

9  50  •  C.  0*0  WILL  S10P  THE  ANALYSIS.') 

9  55  378  CONTINUE 

9  60  READ!  5*  *  )  W1*W2 

965  I  F  <  W  1 .  EG.  0)  GO  10  387 

9  70  GO  TO  160 

975  388  WhI  TEC  6*  390) 

980  390  EORMAK/'  UNITS:  VOLATILITY  MILLIGRAMS/CUBIC  METER'/X 

985  '  HEA1  OE  VAPOhlZATION  KCAL/MOLE ' // ) 

9  90  387  WRI 1  E(  6*  389 ) 

9  95  389  E0RMA1  ('  MORE  DATA?  (  1  =  YES*  0*N0  )') 

1  000  hEADC  5*  *  )  LI 

1005  IE  ( D1  .LI.  0.5)  GO  10  999 

1  010  WRITE  (6*  391) 

1  015  391  EORMAK'  NEW  EILE  NAME  ') 

1  020  READ  C  5*  *  )  AFL 

1  025  CALL  OPEN  <  1*  AFL*  'OUT  PUT  ) 

1  030  GO  TO  1 

1 035  999  END 
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1  050 
i  055 
1  060 
1  065 
1  07  0 
1  07  5 
1  080 
1  085 
1  090 
1  095 
1  100 
1  105  50 
1  1  10 
1115 
1  ISO 
1  1 25 
1  130 


SUBROUTINE  COMP 

COMMON  Y(  100)*  TC  100)#  S4«  S5*A*RjW(4)*V(4)*N1#  IC 

REAL*  8  Y*T#S1#S2*S3#  S4#  S5*  A*  B#  X 

S1  =  0 

S2-0 

S3=  0 

DO  50  K  1=  1#  N  1 

X»- 1.0/CT(Kl)+k( IC)> 

Sl«Sl+X 
S2«S2+X*X 
S3«S3+Y<K1)*X 
CONI  INUE 

B*  (N  l*S3-Sl*S4)/  CN  1*S2-S1*S1) 

A» ( S4-B*S 1 )/N 1 

VC  IC)*N  l*A**2+2.  0*A*B*S  l  +  B*B*S2+S5-2.  0*A*S4-2. 0*B*S3 

RETURN 

END 
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APPENDIX  C 
SAMPLE  RUNS 


hUN 

AN10INE  11:34  10/27/70  TUESDAY  LEA 


UN  11  SI  1  EMPEftAl  Uhk.  CENlILhADE 

FhESSUhE  lOhh 

HOW  MANY  POINTS? 

?  13 

LA1AI  1  Ij  Pli  12*  P2» . 

7  172.1*5.9*  185.3*10.3*  193.4*15.1*  201.0*19.8 


?  211.0*30.1#  218.6*40.4*  227.3*  55.3*  236.7*80.3 


?  251.6*  120.3*  269.3*200.5*  285.0*302.  1*  305.9.502.4*  325.  1*759.2 


DAI  A  TAfcLE?  C  1*=YES1  0«N0  ) 

71 

1 EMPEhATUhfc  PKESSUhE 


1 . 7210D+02 

5.900E+00 

1.8530D+02 

1. 030E+01 

1.9340D+02 

!.510E«-01 

2.  0100D+02 

1  • 980E*  01 

2.  1 100D+02 

3. 010E+01 

2.  I86CD+02 

4 • 040E>  0 1 

2.  2730D+02 

5.530E+01 

2.  3870D+02 

8. 03GE+01 

2.  5160D>02 

1  • 203E*  02 

2.6930D+02 

2. 005E+02 

2.R500D+02 

3.021 E  *  02 

3.  0590D+02 

b. 024E+02 

3.  2510D+ 02 

7. 592Ef 02 

LOG  P-A-E/<C*1> 

A 

7. 06 07 7 D+ 00 

E 

1.89376D+03 

C 

1. 28406E+02 

SI 

ANDABD  DEVIATION 

7. 088  6  b E- 03 

BO 

ILING  POINT 

3.246S0E+02 

29 


i EMPEhATUhE 

PRESSURE 

PERCENT 

EXPERIMENTAL 

CALCULATED 

DIEPERENCE 

172. 100 

5.90000 

5.  73942 

-2. 72167 

185.300 

10.30000 

10.56907 

2.61235 

193.400 

15. 10000 

14.99629 

-0.68683 

201. 000 

19.80000 

20.50014 

3.53603 

21 1.000 

30. 10001 

30.27901 

0.59469 

218.600 

40. 39999 

40. 1 1848 

-0.69681 

227. 300 

55. 30000 

54.55376 

-  1. 34945 

238.700 

80. 30000 

79.82768 

-0. 588 19 

251.600 

120. 30000 

1 19.47359 

-0.68696 

269. 300 

200.50000 

199.09674 

-0.69988 

285.  000 

302. 10010 

301.92554 

-0.  05778 

305.900 

502. 39990 

501.58472 

-0. 16226 

325. 100 

759. 19995 

767.27490 

1. 06361 

TABLE?  GIVE! 

LO  k ES  I #  H I GHES I  AND 

INCREMENT  OP 

TEMPERATURE 

AS  INTEGERS  (0 

#0#0  MEANS  NO) 

?  150#  350#  10 

GIVE  MOLECULAR 

WEIGHT 

?  242.45 

T  EMPERATURE 

VAPOR  PRESSURE 

VOLATILITY 

HEAT  OP  VAPORIZATION 

150.000 

1.81390 

0.  l7E«-05 

20.018 

160.000 

3. 12226 

0.28E+05 

19.546 

170.000 

5. 18220 

0.45E+05 

19.  1  10 

180.000 

8.32321 

0.7 1E+05 

18. 708 

190.000 

12.97603 

0.  1  IE* 06 

18.334 

800. 000 

19.69006 

0. 16E*06 

17.987 

210.000 

29. 15074 

0. 23E^06 

17. 663 

220.000 

42. 19601 

0. 33E+06 

17.36  1 

230. 000 

59.83130 

0.46E+06 

17.078 

240. 000 

83. 24406 

0.63E+06 

16.812 

250. 000 

1 13.8 1445 

0.85E+06 

16.562 

260.000 

153. 12517 

0. 1 IE+07 

16. 327 

270.000 

202.96800 

0.  15E+07 

16. 105 

28  0.000 

265.34741 

0. 19E*07 

15.896 

290.000 

342.48364 

0.  24E>  07 

15.698 

300.000 

436.80786 

0.  30E* 07 

15. 510 

310.000 

550.96069 

0.37E«>07 

15.332 

320.000 

687.78735 

0.  45E+07 

15. 163 

330.000 

850.32446 

0. 55E+07 

15.002 

340. 000 

1041. 79248 

0.66E+07 

14.848 

350.000 

1265. 58472 

0.79E+07 

14. 702 

UNITS!  VOLATILITY  MILL I GRAMS/CUBI C  METER 

HEAT  OP  VAPORIZATION  KCAL/MOLE 


MORE  DATA?  <  1-YES#  O-NO  ) 

STOP 
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UNI  IS*  TEKFERATUhE  CENTIGP.ALE 

PRESSURE  TOF.F 

HOW  MANY  POINTS? 

7  Li. 

DATA!  T1,F1,T2,F2» . 

?  151.6,5.2.  163.0*10.5.  171.0.15.3.  177.0.19.9.  188.  4.32.1 


7  199.  1. 48.  A.  213.3,80.  1,  225. 7,  120.0,  243.  1,202.8 

? 257.8, 302.6,  277.9,501.2,  295.9,755.8 

LAI  A  TABLE?  <  1*  Y  E  S  «  0«N0> 

7  0 


C  LB  55  THAN  1.00000E+02 

1  HE  CONST  AN!  C  IN  LOC  p*A+t/ ( C+T )  IS  ESTIMATED  10  FB 
OUTSILE  THE  RANI  E  100  10  873.  10  pCNTINUE,  SPECIFY  THE 

BNDP0IN15  0 E  ANOTHER  hANLE.  IE  BOTH  ENDPOINTS  APE  1HE 
SAKE,  THE  CALCULATION  WILL  EE  DONE  USINt  THAT  VALUE  EOEi 
C.  0,0  WILL  STOP  THE  ANALYSIS. 

7  50,  100 


LOC  P-A-E/ (C*T  ) 


A  6.21961 L+  00 

P  1.24490D+03 

C  7 • 55994E  *01 

STANDARD  DEVIATION 

FOILING  POINT 

1 EKPEFATUPE 

1. 16  157 B- 02 

2 . 97258  E  +  02 

PRESSURE 

PERCENT 

EXPERIMENTAL 

CALCULATED 

DI E  E  E  FENCE 

151.600 

5. 20000 

5.49919 

5. 75361 

163.000 

10. 500C0 

10.04831 

-4.30181 

17  1.000 

15.30000 

14.83733 

-3.02398 

177.000 

19.89999 

19.55545 

-  1. 73137 

188.400 

32. 10001 

31.92149 

-0. 55611 

199. 100 

*8.39999 

48.72598 

0.67353 

213.300 

80. 10001 

8 1. 37863 

1.59  629 

225.700 

120. 00000 

122.41991 

2.01659 

243. IOC 

202.80000 

205.79420 

1.47643 

257.800 

302.60010 

305.95801 

1. 10968 

277.900 

501. 19995 

498.85498 

-0.46787 

295.900 

755. 19995 

738.93457 

-2. 15378 

TAPLE?  GIVE*  LOWEST, HI  CHEST  AND  INCREMENT  OE  TEMPERATURE 
AS  INTEGERS  (0,0.0  KEANS  NO) 

?  0,0,0 
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DISTRIBUTION  LIST  4 

Agency  Copies 


EDGEWOOD  ARSENAL 

Associate  Technical  Director,  ATTN:  SMUEA-TD-A  1 

Record  Copy,  ATTN:  SMUEA-AGAR  1 

Authors  Copy,  Chemical  Research  Laboratory  2 

Mail  Room,  Research  Laboratories  1 

Director,  Research  Laboratories  1 

Chief,  Medical  Research  Laboratory  1 

ATTN:  SMUEA-RM(3)  7 

Chief,  Veterinary  Medicine  Department  2 

Chief,  Clinical  Research  Department  3 

Chief,  Toxicology  Department  1 

Chief,  Biophysics  Laboratory  1 

Chief,  Industry  Liaison  Office  1 

Chief,  Chemical  Research  Laboratory  1 

Chief,  Analytical  Chemistry  Department  1 

Chief,  Organic  Chemistry  Department  1 

Chief,  Physical  Chemistry  Department  1 

Chief,  Physical  Research  Laboratory  1 

Chief,  Dissemination  Research  Department  1 

Technical  Support  Directorate 

Chief,  Field  Evaluation  Division  1 

Technical  Information  Division 

ATTN:  SMUEA-TSTI-TP  2 

ATTN :  SMUEA-TSTI-TR  1 

ATTN:  SMUEA-TSTI-TL  10 

Director,  DDEL 

ATTN:  SMUEA-DPP  1 

ATTN:  SMUEA-DDW  1 

Director,  WDEL 

ATTN:  SMUEA-W  1 

ATTN:  SMUEA-WCP  1 

Systems  Analysis  Office,  ATTN:  SMUEA-OS  1 

Patent  Advisor,  Office  of  Counsel  1 

Foreign  Intelligence  Office,  ATTN:  SMUEA-POF  3 

Director,  USAMUCOM,  Operations  Research  Group  1 

Commanding  Officer,  US  Army  Technical  Escort  Center  1 

CO,  Quality  Assurance  Directorate,  ATTN:  SMUEA-QAE  1 

DEPARTMENT  OF  DEFENSE 

Defense  Documentation  Center  20 

Cameron  Station 
Alexandria,  Virginia  22314 
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